Short and long range correlated motion observed in 
colloidal glasses and liquids 

Eric R. Weeks W, John C. Crocker ^ and D. A. Weitz^ 3 ) 

^Physics Department, Emory University, Atlanta, GA 30322 
E-mail: weeks@physics.emory.edu 

( 2 ) Department of Chemical and Biomolecular Engineering, University of 
Pennsylvania, Philadelphia, PA 19102 

(•^Department of Physics and DEAS, Harvard University, Cambridge, MA 02138 

Abstract. We use a confocal microscope to examine the motion of individual 
particles in a dense colloidal suspension. Close to the glass transition, particle motion 
is strongly spatially correlated. The correlations decay exponentially with particle 
separation, yielding a dynamic length scale of 0(2 — 3d) (in terms of particle diameter 
a). This length scale grows modestly as the glass transition is approached. Further, the 
correlated motion exhibits a strong spatial dependence on the pair correlation function 
g(r). Motion within glassy samples is weakly correlated, but with a larger spatial scale 
for this correlation. 



PACS numbers: PACS: 61.43.Fs (glasses), 64.70.Pf (glass transitions), 82.70.Dd 
(colloids), 61.20.Ne (structure of simple liquids) 
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1. Introduction 

The viscosity of a glass-forming material increases rapidly by many orders of magnitude 
as it is cooled, without any corresponding structural change to account for the viscosity 
change [U |2]. Adam and Gibbs suggested that a growing dynamic length scale may 
relate to the viscosity growth [3 . Their idea has been interpreted in several ways 
[H IH El HI El but no experiment has been able to directly observe dynamic 
length scales, nor is it clear what form these length scales would take. Indirect 
evidence is provided by experiments which locally perturb materials in a variety of 
ways, and observe the relaxation response of the material and its dependence on the 
scale of the perturbation j^]. Further evidence for a dynamic length scale comes from 
experiments done in thin films or small pores [TO] , but these experiments are unable to 
directly observe the nature of any correlated motion. Recent simulations have examined 
correlation functions in systems of Lennard- Jones particles [17] and hard spheres 
finding evidence for spatial correlations and a possible dynamical correlation length scale 

[13 IUJ. 

We study a system of colloidal particles which interact only via repulsive forces, 
and which have a glass transition as their concentration is increased. We use a confocal 
microscope to track the motions of several thousand particles for several hours, long 
enough for most particles to undergo nontrivial displacements. The mobilities of 
these particles (the magnitudes of their displacements) are correlated over distances 
of ~ 2 — 3a, in terms of the particle diameter o~. These correlations reflect large- 
scale cooperative rearrangements of particles seen previously [TTJ H2J E3 UTJ I2H|- We 
also examine the correlations between the directions of particle displacements. While 
the directions are correlated, this does not appear to be as strong a signature of the 
rearrangements. These measurements are direct experimental studies of the nature of 
the long-range correlations present near a glass transition, and indicate that rearranging 
regions of particles are composed of particles with large displacements, but which do 
not all move in similar directions; rather, these regions are internally rearranging [28j . 

2. Experimental Methods 

Our system consists of a suspension of colloidal poly-methyl-methacrylate, sterically 
stabilized and dyed with a rhodamine dye [T^j. These particles are slightly charged, 
and have a hard-sphere diameter a = 2.36 /im and a polydispersity of 5%. They 
are suspended in a mixture of cycloheptylbromide and decalin which nearly index 
and density matches the particles. The samples are prepared with a constant volume 
fraction <ft and sealed into microscope chambers. We observe crystallization occurring at 
cf) = 0.42, slightly lower than the value expected for hard spheres (4>hs — 0.494) [2HI22]- 
Samples with > <p g « 0.58 do not form crystals within the bulk even after sitting at 
rest for several months, allowing us to identify <p g as the glass transition volume fraction, 
in agreement with previous work |2l| I25j . Prior to observation, samples with crystals in 
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them are shear-melted by means of a stir bar; there is a reasonable separation of time 
scales between the decay of transient flows caused by the stirring (< 20 min) and the 
onset of crystallization for these samples (> 5 - 10 hours, as determined by methods 
described in j"*"7j). 

We use confocal microscopy to observe the particle motion. A single three- 
dimensional image is acquired in 10 s, and over the course of an experiment we acquire 
several hundred images. As the mean square displacement curves show in figure EJa), 
particles do not move significant distances on a time scale of 10 s. We post-process 
the data to determine particle positions with an accuracy of 0.03 /im horizontally 
and 0.05 /im vertically [13]. Each three-dimensional image is 69 x 65 x 14 /zm and 
contains several thousand particles. We focus at least 25 /mi from the coverslip to avoid 
interference from the wall, and in fact the observed motion appears isotropic. 

3. Results 

To characterize the behaviour of the samples, we calculate the particles' mean square 
displacement (Ax 2 ), shown in figure IHa) for five samples. These curves all show a 
plateau, due to cage-trapping: each particle is confined in a cage formed by its neighbors. 
At longer times, (Ax 2 ) shows an upturn, indicating that at least some particles have 
moved. For the supercooled fluids (thick lines, < (f> g ), previous work showed that 
these motions correspond to a small subset of particles undergoing cage rearrangements 
|11[ IT2"1 ITT| 123*] . These particles move significantly farther than the majority of the 
particles, and thus distributions of the particle displacements show broad tails at these 
time scales [TTJ "T21 [J7J [23] • This is quantifiable by the nongaussian parameter 

°^ = 3^ - 1 W 
where the moments of Ax are calculated from the measured distributions of the 
one-dimensional displacements Ax. a.2 is zero for a gaussian, and larger when the 
distributions are broader than a gaussian. We plot «2 i n figure H*[b), finding that for 
supercooled fluids (thick lines) a% has a peak corresponding to the end of the plateau 
of (Ax 2 ), due to the presence of the anomalously mobile particles [TT]. For colloidal 
glasses [thin lines in figure da)], the upturn in (Ax 2 ) is due to aging, and occurs at a 
time scale At which varies with the time since sample preparation *2H]- It is unclear if 
the motions responsible for the upturn are due to cage rearrangements [TT] . 

Intriguingly, in supercooled colloidal fluids, the motion of these cage-rearranging 
particles is spatially localized, and the peak of a 2 (At) corresponds to the existence of 
large clusters of these particles all moving simultaneously [TT]. In fact, the positions of 
mobile particles appear in localized clusters over a range of time scales At, and the At 
dependence of the typical cluster size appears qualitatively similar to the At dependence 
of a 2 [UJ. 

We wish to determine how these previously observed clusters of cage-rearranging 
particles are manifested in correlations of the motion of individual particles. In a 
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Figure 1. (a) Mean square displacement for colloidal "supercooled" fluids (thick lines) 
and colloidal glasses, with volume fractions <j> as labeled, (b) Nongaussian parameter 
«2 for the samples shown in (a). (c,d) Vector correlation function Ss(At) and scalar 
correlation function S$ u (At) for the samples shown in (a), using Ar corresponding to 
the first peak of the peak of g(r) for each sample. 



dilute suspension of particles, particles have only hydrodynamic interactions, and the 
correlations between the displacement vectors u of two particles should decay as 1/Ar 
with increasing separation Ar [TBI EE 120] • This result would still hold true for dilute 
tracers in any homogeneous viscoelastic medium [3UJ ETJ. However, it seems unlikely 
that our samples are homogeneous on the length scales we consider; our imaging window 
has a width of only ~ 25 particle diameters. The direct interparticle forces (their electric 
charge, and steric repulsion) are likely to be as important as hydrodynamic interactions. 

Thus, we do not necessarily expect 1/Ar decay of correlations, especially given 
the apparently cooperative and localized nature of the cage rearrangement motion. 
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For example, a simulation of hard spheres at large volume fractions found exponential 
decay This suggests particles may rearrange by translating together in a group. 
One possibility for this group motion is that cage rearranging particles move in parallel 
directions fH^JEHEl- ^ n alternate possibility is that mobility may be correlated over 
long distances. Mobility is the magnitude of particle displacements u — \u\, sometimes 
with the average subtracted off (5u = u — (u)). In fact, the previous observations 
of clusters of anomalously mobile particles are a direct indication that mobility is 
correlated; simulations have provided evidence that there is a correlation length scale 
• Correlations of mobility are consistent with both collectively translating regions, 
and also internal rearrangements within localized regions. 

To study these possibilities, we compute two correlation functions from our data, 
one using u and one using the mobility 5u [TBI E01 \U\ • We define 



where the time scale At is used to define the displacements u and mobility 5u. For 
both formulas the numerator average is over all pairs of particles i, j with separations 
Ar, and the denominator average is over all particles fZT\ 122] . The normalization of 
both of these functions is chosen so that perfectly correlated motion {ui = Uj for all 
particles corresponds to Sg = Ss u = 1, perfectly anticorrelated motion (vii = —Uj) 
corresponds to = —1, and uncorrelated motion gives S$ = Ss u = 0. Note that 
Ssu measures fluctuations of mobility, so either anomalously mobile or anomalously 
immobile particles would show positive correlation. By examining the Ar dependence 
of these functions we learn about spatially correlated motion, and by examining the At 
dependence we study the correlated motion seen at different time scales. 

To study the At dependence of S$ and S$ u , for convenience we choose a fixed value 
Ar = Ar n „, for each data set, set by the first maximum of the pair correlation function 
g(r) (which is slightly larger than at 2a, and depends slightly on <ft, due to the charging 
of the particles). Pairs of particles with separation Ar nn are nearest neighbors and are 
strongly correlated due to their close proximity, as will be shown later. 

We plot Sa(At,Ar nn ) in figure dc) and Ss u (At,Ar nn ) in figure UJ^d) . For the 
supercooled fluids (thick lines), the behaviour of these two functions is dramatically 
different. is nearly constant as a function of At, whereas Ss u is small at short lag 
times, rises to a maximum, and then has a downturn. Strikingly, the At dependence 
of Ssu is similar to that of the nongaussian parameter a 2 , which can be seen for three 
supercooled fluid samples by comparing the thick lines in figure Hfb) and figure Efd). 
While Ss u reaches its maximum at slightly larger values of At, it appears that the 
anomalously mobile particles (which cause the increase of a.2) are directly responsible 
for the increasing correlation measured by Sg u - While the time dependence of the 
vectorial correlations S$ is much weaker than for the mobility, there is a slight rise at 
larger At seen in figure Hfc) [IB]. It has been conjectured that S^(At — > 00) should be 




S Su (Ar, At) 




(2) 



(3) 
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nonzero, as short-lag-time interparticle correlations will always provide a contribution 
to (ui ■ Uj) even if particles are uncorrected at longer lag times [T^]. This may explain 
why we do not see a downturn in S$ at large At. 

Figures [Tfc, d) also show that at short lag times, > Ss u , although at the cage- 
breaking time scales, the opposite is true. These results, along with the similar lag time 
dependence of a 2 and S$ u , suggest that the cage rearrangements are due to mobility 
correlations, rather than a strong directional correlation. The implied picture is that 
regions of cage rearrangements are composed of highly mobile particles, which move in 
many directions. While the motions of neighboring particles are somewhat directionally 
correlated, the cage rearrangements reflect regions of internal restructuring rather than 
large-scale cooperative translations. 

Further confirmation of this picture comes from the volume fraction dependence 
of the correlation functions, seen in figure ^c,d). As increases towards the glass 
transition at (fi g ~ 0.58, changes only slightly, while Sg u changes dramatically, again 
similar to the behaviour of a.i- The growth of Ssu (and non-growth of Sn) indicate 
that as the glass transition is approached, mobility correlations become increasingly 
important, while the directional correlations remain virtually unchanged. 

The results are harder to interpret for glassy samples (thin lines in figure due 
to the aging of the system. In an aging system, the properties of the samples depend 
on the time since the sample was prepared [2H]; for example, the increase in (5x 2 ) at 
large At [figure ^a)] moves to longer At as the sample age increases. (The sample age 
is defined as the time since stirring the sample prior to the start of the experiment. To 
minimize the effects of aging for these data, data acquisition was started 5-10 hours 
after stirring the sample, which reliably initiates the aging [2*3 ■) As seen in figure 
IHc,d), any cf> dependence for the correlated behaviour in glassy samples is unclear. The 
correlation functions are nonzero due to correlated motion of particles within cages; 
however, the increases seen in and Ss u at large At are probably due to the slight 
motions responsible for aging [29 . 

To examine the spatial character of the highly correlated particles occurring at the 
cage-rearrangement time scales, we study the Ar dependence of these functions fixing 
At to maximize Sg u [figure EJd)]. The correlation functions are shown for a supercooled 
fluid (<f) = 0.56) in figure 12(a), indicated by the thick lines. Both functions oscillate, with 
especially strong oscillations seen in Su (thick solid line). Strikingly, these oscillations 
coincide with oscillations of the pair correlation function g(r), shown in figure 12(c) . For 
example, pairs of particles with separations corresponding to the first peak in g(r) are 
nearest neighbors, and are the most strongly correlated. This agrees with our earlier 
work which found that the pair correlation function g(r) influences the motion of pairs 
of particles 28] . In fact, comparing Ss u and Sn shows that the correlation of the mobility 
of pairs of particles is less sensitive to g(r), whereas the correlation of their directions 
is more sensitive to g(r). This is further evidence that cage rearrangements involve 
regions of particles with high mobility, although within those regions the directions of 
the displacements of individual particles are not as strongly correlated, but rather are 
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Figure 2. (a) Correlation functions for <j> = 0.56 (a supercooled fluid close to the glass 
transition), with At = 2000 s. (b) Correlation functions for <f> = 0.62 (a glass), with 
At = 10000 s. (c) g(r) for the data shown in (a); g(r) is similar for the data shown in 
(b). 



influenced by g(r). Both simulations of Lennard- Jones particles and of hard spheres 
saw similar oscillations in spatial correlation functions, but the relationship with g{r) 
was unclear as the systems had higher polydispersity than our samples [T7| ITol 122]. 

To further examine the correlations in the directions of motion of pairs of particles, 
we decompose into longitudinal correlations (correlations of the displacement along 
the direction of the separation vector Ar*jj) and transverse correlations (correlations of 
the component of the displacement vector perpendicular to Ar^). Two new functions 
are defined as: 



Su L (Ar, At) 
SW(Ar, At) 



(ff • iff) 



(4) 

(5) 
The 

Sul + SsT- Rij is a unit vector pointing from 



T 



U, 



— ufRij, and uj 



u ; 



uf Rij- 



where uf 

normalization is chosen so that Sq 
particle i to particle j, and the dot products for the u L terms are taken so that two 
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Figure 3. Semilog and log- log plots of S$ (solid lines) and Ss u (dashed lines) at two 
different At as indicated, for <f> = 0.56, a supercooled fluid. The dotted 1/Ar lines are 
drawn as a guide to the eye in the log-log plots. The dotted fit lines shown in the lower 
left plot have decay lengths £su — 8.5 /im and £3 = 6.1 fim. The downturn at long At 
is not fitted, and is described in the text. 

vectors pointed in the same direction correlate positively. These two functions are also 
plotted in figure Efa) (thin lines). The oscillations of Sg are almost entirely due to 
the contribution from S$l (thin solid line), indicating that the longitudinal motion is 
most sensitive to g(r). In fact, the longitudinal correlations of the motions of nearest 
neighbor particle pairs are suggestive of "string-like" motion seen in simulations 
However, S$i also shows a slight anticorrelation at Ar pa 3.2 /im, corresponding to the 
first minimum of g(r). This indicates that pairs of particles separated by ~ 1.5a are 
actually more likely to move in antiparallel directions (towards or away from each other), 
despite their strongly correlated mobility. The transverse component S^t (thin dotted 
line) has a slight increase closer to the first minimum of g(r), but otherwise shows little 
dependence on g(r). 

At larger separations, particle motion is still less correlated. If the particles were 
in a dilute suspension and long range interactions were solely due to the hydrodynamic 
behaviour of the solvent, the correlation functions should decrease as 1/ Ar. This would 
also be the case for a homogeneous viscoelastic medium j20j, although the localized 
rearrangements observed previously shows that motion in these dense colloidal samples 
is spatially inhomogeneous fTJ To check this, we plot the correlation functions 
on semilog and log-log axes in figure El for = 0.56 at two different time scales. The 
top right graph shows that the decay is close to 1/Ar for short time scales; we find 
similar behaviour for all liquid samples at time scales At < 100 s. The behaviour 
characteristic of longer time scales relevant for cage rearrangements is shown in the 
bottom graphs. On the semilog graph, the decays follow straight lines for both functions 
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Figure 4. Volume fraction cj> dependence of the dynamic length scales, for colloidal 
liquids and glasses. The length scales are determined by a fit to the exponential decay of 
correlation functions such as those shown in figure|3| the error bars indicate variability 
in determining the length scales at different Ai's. (The data points for 4> = 0.52 and 
4> = 0.56 are horizontally offset from each other slightly for clarity.) 



out to Ar « 25 /im, indicating exponential decay with characteristic length scales 
£,s u = 8.5 /im = 3.6a and = 6.1 /im = 2.6a. For this sample, the decay lengths 
vary in the range £$ u = 7.6 ± 0.9 /im and = 6.7 ± 0.9 /im when At is varied from 
50 - 5000 s. The exponential decay demonstrates that for dense colloidal suspensions 
on these length and time scales are not behaving as continuum viscoelastic materials 



At large Ar, a downturn is seen in the correlation functions (bottom graphs in figure 
0J). This is due to a counterflow, similar to what was seen in simulations [TH| I26j. This 
counterflow has been previously interpreted as the medium's response to the transient 
motion of a particle, although it is less clear how this may apply in the case of localized 
motions [TSJI^n]. The counterflow may be due to the presence of the coverslip (>30 /im 
away), although the correlation functions do not change significantly when the data 
is split into two subsets, one closer to the coverslip and one farther. The counterflow 
results in a slight anticorrelation in both and S$ u at Ar « 40 /im (not shown). In 
all cases the exponential decay occurs over a decade in Ar before the counterflow cuts 
off the correlation. At still larger length scales ( Ar ^> £g u ) , the correlation functions 
must decay as 1/Ar, as on a sufficiently large length scale the colloidal suspension will 
appear homogeneous. However, we do not have data at large enough Ar to see 1/Ar 
decay, and at large separations the amplitude of the correlation functions would be quite 
small and difficult to measure. 

To look for a possible growing length scale, we extract the correlation decay lengths 
for Sn and Ss u in samples with different volume fractions, with results shown in figure 0] 
The decay lengths increase only slightly as (p g is approached. Within the error bars, £g 
is consistent with a constant value £g ~ 3a. £$u doubles over the range = 0.46 — 0.56, 



EH EI]. 
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and this is sufficient to account for the increase in the size of clusters of mobile particles 
seen seen in previous work [TI]; however, we see no evidence in our data for or against 
a divergence of £$u at <$>g- Intriguingly, this slight increase in £su gives new insight into 
the 4> dependence of S$ u for the supercooled fluids, seen in figure Efd). If we assume 
Ssu(Ar) = A^At) exp(— Ar/^su), we find that the maximum value (with respect to 
At) of A^At) is a constant, approximately 0.75 ± 0.05, independent of 0. Thus, 
the increasing height of Ss u as increases, seen in figure W[<i), is primarily due to 
the increase in £g u and thus an increasing value of exp(— Ar nn using the nearest 
neighbor separation Ar nn . 

We can compare our length scales and with length scales obtained from 
computer simulations of glass-forming systems. The values we find for = 0.56 
(&z ~ £su ~ 3cr) are slightly larger than that seen in simulations of hard spheres 

= 2.3cr, £su = the simulations had a higher polydispersity, and thus may 

have been farther from the glass transition at = 0.56 Simulations of Lennard- 
Jones particles did not notice exponential decay, and instead computed a length scale 
by calculating (in our notation) 



The term in parenthesis is qualitatively similar to the nongaussian parameter a-z, and 
over the range of liquids shown in figure 0] varies from 0.35 (at = 0.46) to 0.70 (at 
= 0.56). We calculate £' for our data by using the form A^(At) exp(— Ar/£su) for Sg u 
with our measured values for A and £s u , and using the value of At which maximizes 
Ssu (and thus maximizes £'). We find that £' increases from 0.9 - 4.1 /im « 0.4<r — 1.7cr 
as the glass transition is approached. These values are larger than those seen in the 
simulation (0.05 - 0.32 in Lennard- Jones units) |23j . 

Earlier work has seen evidence that structural properties of the sample are slightly 
correlated with particle mobility [2HIE21- This then suggests that the exponential length 
scales seen in the correlation functions may relate to spatial correlations of structural 
properties. We have computed a spatial correlation function Ssv similar to S$ u (3), 
except looking at fluctuations 5V of the Voronoi volume V of each particle. The Voronoi 
volume is a geometric way of partitioning space, so that each particle claims the volume 
that is closer to its center than to the center of any other particle. Particles with slightly 
larger Voronoi volumes are, in some sense, seeing a slightly smaller local volume fraction, 
and this is correlated with a slightly enhanced mobility [28J. The correlations Ssv do 
decay exponentially (not shown), with a length scale of 5.7 — 6.3 /im for the four liquid 
samples we consider (~ 2.5<r). Strikingly, we see no clear volume fraction dependence 
of this length scale. Furthermore, this length scale is quite similar to the length scales 
seen for S^. The lack of a clear connection to the length scale for mobility (£su) suggests 
that the correlation between the structure and the mobility, while present, is not the 
entire story; the cooperative motions are influenced by more factors than just the local 
volume [USES]. 




(6) 



Short and long range correlated motion observed in colloidal glasses and liquids 11 

For the glassy samples, the behaviour is strikingly different. As discussed before, 
our glassy samples are aging, although for the duration of the experiment, little aging 
occurs. (We consider samples with a large age t w = 3-10 4 s and time scales At <t w .) For 
our data, and for lag times At < 5000 s, almost all particle motion consists of particles 
moving back and forth within their cages. and Ss u have a low amplitude, as shown 
in figure ^c,d) and figure I2fb), indicating that most of the motion is uncorrelated. For 
the glassy samples, S& > Ss u , similar to the liquids at small At. Figure Efb) also shows 
that the correlation functions oscillate much less, and thus the local structure [g{r)} 
has a only a minor role for glasses. Another intriguing difference is that S^t > S^l, 
indicating fluctuations in the transverse directions are more strongly correlated than in 
the longitudinal direction. 

The correlated motion that is present in the glasses, while of small amplitude, 
is rather long-ranged. Similar to the supercooled fluids, the correlation functions for 
glasses show exponential decay over a range of time scales At. The length scales for 
the decay are plotted in figure EJ They are significantly larger than the length scales 
associated with the supercooled fluids; for the glasses, £„,<5« = 8 — 20 fim « 3.5a — 8.5a. 
The data shown in figure |U also indicate that the scalar correlation lengths £$ u are 
noticeably larger than the vector correlation lengths a trend which is the opposite 
that of the liquids. The existence of this long-range correlation is sensible, given that 
the glasses are more densely packed: while most motion is localized and uncorrelated, 
particles can move slightly if there is long-range cooperation, even if they then move 
back (in order that there be no rearrangements). At present, however, it is unclear if 
these correlation length scales correspond to the intrinsic motion of the frozen glassy 
system, or if they are connected with the aging process [29 . The previous observations 
that clusters of "mobile" particles are smaller is consistent with the low amplitude of 
the correlations [TTJ . 

Unlike the correlated motion in glasses, the structural correlations in the glassy 
samples are not long-ranged. The length scale for Sgv, the correlation function of the 
fluctuations of Voronoi volumes, stays comparable to that of the liquids (£sv ~ 5.9 — 7.7 
(Ltm for the glasses). While there is very little change of this length scale with 0, we 
do note that the two samples with the largest values of £sv correspond to the same 
two samples with the largest values of again suggesting a connection between the 
structure and correlated particle motion [3"^| . 

4. Conclusions 

By measuring the positions of several thousand particles over several hundred time 
steps, we have found that particle motion is correlated over distances of 3-4 particle 
diameters. In particular, particles undergoing cage rearrangements have anomalously 
large displacements, and these are highly spatially correlated. We note that the largest 
degree of correlation exists at the time scales corresponding to cage rearrangements, 
but that our observations do not preclude the sample from acting differently at very 
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long time scales. All of our data is taken at time scales At < r a , so we are not able to 
quantify the behavior on time scales corresponding to alpha relaxation. 

At time scales similar to the cage rearrangement time scale, the correlations decay 
exponentially with separation Ar. At very large Ar, we would expect the correlations 
to recover a 1/Ar dependence as predicted for homogeneous viscoelastic media [HOI EI], 
at any time scale. This dependence is seen at short time scales in our data (figure 
EJ), suggesting that the instantaneous response may be more continuum-like [33]. This 
is sensible; in such short time scales, At <C r a , and so these motions are not related 
to dynamical heterogeneities, and do not lead to flow of the sample. Rather, these 
are simple small amplitude affine deformations. It is intriguing that on the cage- 
rearrangement time scales, our image volume is not large enough to see the 1/Ar 
behavior, suggesting that the appropriate coarse-graining length scale is bigger than 
~ 25(7. 

The differing behaviour of the two correlation functions we examine suggests that 
these rearrangements are composed of regions of mobile particles (particles with large 
displacements) , but that the motions involve particles moving in many directions, rather 
than all of the particles moving in one coherent direction. This is especially clear as 
the correlations between the directions of motion of particles (Sy strongly depend on 
the pair correlation function; within a group of mobile particles, their motions are often 
strongly directionally correlated, but also frequently not directionally correlated. This is 
consistent with particles moving in necklace-like loops, for example, as seen in Ref. |16j . 
and the mixing motions described in Ref. [28J. 

The correlations become increasingly long-ranged as the glass transition is 
approached, as seen in two new correlation lengths £su and This suggests that 
rearrangements involving regions consisting of a small number of particles become 
difficult or impossible as the volume fraction increases, which would explain the growth 
of the viscosity as the glass transition is approached. The size of these regions, quantified 
by the two correlation lengths £g and £s u , are direct experimental evidence for dynamical 
length scales near the glass transition. 
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